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ABSTRACT 


Ar experimental and theoretical analysis of a rotating system with 
flexibly supported, viscousiy damped bearings is described. The theo- 
retical develcpment is based on the differential equations of motion of 
the mass of the rotor and the two shaft bearings. Tests were conducted 
and compared with caiculated response of the system. Results indicate 
fair agreement The discrepancies are attributed to limited accuracy 


of test measurements and/or erratic operation of the viscous damping 
dashpots. 
This work was conducted at the U. S. Naval Postgraduate School, 


Monterey, California. 
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1 INTRODUCTION 


Tre vroblem cf vibration in rotating systems has commanded the 
attention of theoreticians and practical investigators for many years. 
With the advent of high speed marine turbines and turbo-jet aircraft 


engines renewed interest has been shown in this field. 


it is weli krown that at certain speeds of rotation small imbalances 
tend to set up vibrations of large amplitudes in a rotating system. 
Probl, who has made several contributions in this area, а , wrote, 

The existing methods for determining critical speeds are 

subject to the following limitations. On one hand the methods 

that are general, i.e., that permit the calculations of higher 

speeds as well as the fundamental, involve calculations so 

complicated as to be impractical for any but the simplest of 
rotors. On the other hand, the methods for which the calculations 
are comparatively simple, such as the more familiar methods of 

Rayieigh and Stedola, lack generality in that critical speeds 

other than the fundamental cannot be definitely determined. 

This criticism may have been well taken at the time. The large 
number of prime variables in any but the simplest system tend to confine 
most investigations to a very narrow scope; for example, "On Critical 
Speeds of a Shaft Supported by a Ball Bearing", by T. Yamamoto, which 
appeared in the Journal of Applied Mechanics of June, 1959. Doctor 


Yamamoto investigated two kinds of critical speeds induced by a slight 


difference in the diameter of the balls in their races. 


A growing appreciation of modern computer capabilities has led to 
the development of increasingly general treatments such as, "Natural 
Frequencies of Nonuniform Beams on Multiple Elastic Supports", by 
К. A. DiTaranto cf the Defense Electronics Products Division of Radio 
Corporation of America, (31. An especially general and pertinent 
method ws described recently by Messrs. Koenig, Guenther and Love joy 
of the Engineering Analysis Section, Allis-Chalmers Manufacturing 
Company, 14]. A concluding, third part of this method, as yet unpub- 
lished, is to include experimental results and correlation with the 


analytıcal method described 
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The study presented here investigates the transverse motion of 
a rotor оп а flexible shaft. The motion of the two flexibly supported 


end bearings, which are subject to damping, is also studied. 


The study includes both a practical investigation conducted with 
a machine designed specifically for this purpose and a theoretical 


analysis developed along parallel lines. 


The theoretical analysis begins with the basic equations of motion. 
This approach was considered preferable to a mechanical application of 


one of the existing schemes of analysis. 


The analytical model was constructed with the assumption that the 
rotor was the only rotating member. The shaft was considered to be 
massless and torsionally rigid. The speed of rotation was taken to 
be constant. Damping was assumed to be a linear function of the 


velocity. 


This work was conducted at the U. S. Naval Postgraduate School, 


Monterey, California during the 1961-1962 academic year. 


The author wishes to express his appreciation to Professors 
T. H. Gawain and M. H. Vavra and to the supervisory and laboratory 
personnel of the Department of Aeronautics of the U. S. Naval Post- 


graduate School. 
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2 TABLE OF SYMEOLS 


general constant 
Fortran array 


rotor position measured from left-hand support, in. 


(1, 2, etc.) Fortran physical constant 
(A, E) Fortran damping constant 


symmetrical damping constant, lbs /(in./sec.) 
displacement, in 
mean spring diameter, in. 
wire diameter, in. 
force, ibs. 
shear modulus, ibs./sq. in. 
gravitational constant, 386 in./sec. 
dismetral moment of inertia, in. 
i 
imaginary, (-1)^* 
spring constant, 155. /1п. 
shaft iength, in. 
mass. о a. 
number of working coils of spring 


first critical speed, RPM 


load, ibs. 

general constant 

general constant 

time, sec. 

weight, lbs. 

rim weight, Ibs. 

displacement of rotor c.g. from polar axis, in. 
displacement, in. 

velocity, in./sec 

acceleration, Тсс. 


general constant 





Table of Symbols 


5 complex root 
Е rover eccentricity, in. 
А damping factor, - 
B mass ratio 
+ 

un Py poz 
5 ы Hz 

3,1416 
@ phase angle, degrees 

angular frequency, radians/sec. 

Subscripts ° 
1 left-hand bearing complex 
ZR rotor 
2 right-hand bearing complex 
d damping 
I inertia 
T total 
Abbreviations: 

c.g center of gravity 
cps cycles per second 
in. inches 
lbs. pounds 
шп. mass units 
RPM revolutions per minute 
sec. seconds 
sq.in. square inches 





3. THE PHYSICAL SYSTEM 


A schematic drawing of the physical system is shown in Fig. 1. 
Detailed drawings from which the critical speed test rig was constructed 


are included in Appendix A. 


3.1 Critical Speed Rig Design Concept 


Bearings and restraints of high-speed rotating systems should be 
designed with great care. Bearings are normally confined as rigidly 
as possible in order to reduce the degrees of freedom of the system and 
to eliminate critical speeds below the operating speed range. However, 
shaft and bearing flexibility are inevitable due to the elastic nature 
of materials, oil film 'stiffnesses", necessary clearances and machining 
tolerances. Similarly, imbalances due to imperfections and lack of 
homogeneity of shaft and rotor materials are always present to some 
degree. Even after exhaustive balancing procedures have been carried out 
some imbalances remain and contribute harmonic excitations in rotating 


systems. 


The term impedance is sometimes used in describing the character- 
istics of spring and damping device combinations which act at the supports 
of a system. As noted by Guenther and Lovejoy, 177 Support mechanical 
impedance can be, and is often, expressed as a function of not only the 


static properties of a system but of the operating speed as well. 


In the case where support is provided by a simple spring, support 
stiffness is merely the familiar spring constant (force per unit 
deflection). The ideal mechanical spring acts as a conservative member 
in the system; energy stored in compression or tension is returned to 
the system when the spring returns to its unstrained state. Hysteresis 
losses present in actual springs are generally small and can be disre- 


garded. 


For the simple, damped, spring-mass system as discussed in most 
basic vibration texts, e.g., os a dashpot is shown as a device which 
develops forces opposing motion, the force so developed being propor- 
tional to the velocity. A dashpot acts as a non-conservative member 


dissipating energy from the system. 
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Both the spring and the dashpot influence the critical speed(s). 
The dashpot governs the amplitudes associated with critical speed 


vibrations 


The critical speed test rig was designed to provide a vehicle for 
practical investigation of the effects of support flexibility and 
viscous damping on the critical speeds and mass vibration amplitudes of 


a rotating system. 


The salient points of the design with respect to its particular 
service are: (a) compressed air drive as opposed to mechanical drive 
which would constrain the shaft during vibration, and, (b) ability to 


control dashpot stiffness by means of a suitable valve. 


Choice in arrangement of the three support frames, the several 
rotor and spring sizes available, range in dashpot stiffness and choice 
of rotor position on the shaft afford a broad selection of configurations 


possible of investigation with this test rig. 


3.2 Test Configuration and Major Components 
3.2.1 Drive Disk and the Bearing and Linkage Complex 


Drawings of these components are shown in Appendix A. Pertinent 
properties are listed in Table I, Physical Properties of Components. 
Figure 2 shows the drive assembly and Fig. 3, the bearing and linkage 


complex. 


и оз 
Figure 3 shows springs mounted in support of one of the shaft 
bearings. Details of the various springs are shown in Appendix A and 


listed in Table I. 


Spring constants were calculated by formula, (el. 


P = load, ibs. 
F = deflection caused by given load, in. 
G = shear modulus, lbs./sq. in. 
d = wire diameter, in. 
D = mean diameter of spring, in. 
N = number of effective working coils 
А 

== => 115. /1п. 
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Light Spring 


11.4 x 10° x 2.07 x 104 


с Та 





= Zn 


Intermediate Spring 


-4 

Spring constants were also determined experimentally. The spring 
and its hoider were set vertically over a slot in a fixed surface. A 
small wire attached to the top holder was led along the axis of the 
spring and through the slot in the table. The unloaded length of the 
spring was measured with a vernier caliper. Weights of known sizes 
were attached to the free end of the wire. The length of the spring 
was measured for each increment of step-wise loading, and similarly 
during unloading. Results are listed in Table II, Calibration of Light 
and Intermediate Springs. Results are shown graphically in Figs. 4 and 5 
for light and intermediate spring respectively. Calculated and measured 
values of spring constants are listed in Table III, Measured and Calcu- 


lated Spring Constants, for purposes of comparison. 


The lack of agreement in the calculated and test values for the 
light spring indicates that friction in the holders may have had a strong 
influence on the experimental results. Further testing of the spring is 
advisable. Both Fig. 4 and 5 show pronounced hysteresis losses are 
present during compression and relaxation of the springs. The slopes of 
the figures, from which the response was calculated, were not strongly 
effected by hysteresis. The experimental value of the intermediate spring 
constant, in good agreement with the calculated value, was used in all 
subsequent calculations Only the intermediate spring was used in tests 


for critical speeds. 
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Table II 


CALIBRATION OF LIGHT AND INTERMEDIATE SPRINGS 


LIGHT SPRING INTERMEDIATE SPRINGS 
Wt.(lbs.) Deflection(in.) Wt.(lbs.) Deflection(in.) 
0.0 0.0 0.0 0.0 
1.0 0.053 KO 0.061 
2.0 0.062 2:50) 0.109 
330 0. 066 3.0 Ори Ба 
4.0 0.078 4.0 0.238 
5.0 0.083 5.0 0.290 
6.0 0.091 6.0 0. 382 
7.0 goo 70 0:513 
8.0 0.135 8.0 0.597 
9.0 0.146 7.0 0.550 
10.0 0.168 6.0 0.506 
9.0 0.149 5-0 0.447 
8.0 0.143 4.0 Û0. 365 
7.0 0 153 3.0 0.281 
6.0 07122 220 0.214 
5.0 0.110 1.0 0.142 
4.0 0.101 0.0 0.101 
Boo 0.087 
2-0 0.082 
1.0 0.070 
0.0 0.062 
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Table III 
MEASURED AND CALCULATED SPRING CONSTANTS 


ж 
Spring Measured Calculated % Variation 
Constant Constant 
Light 14.5 17.5 2075 
Intermediate 105.5 197.0 ПА 


* 
based on measured value 


3.2.3 Shaft 
Dimensions and physical properties are shown in Appendix A and 


Table I. 


The shaft was checked for set or eccentricity by rolling on a 
precision flat surface. Variations noted were so slight that they did 


not warrant correction. 


Static shaft spring properties were calculated for midspan loading 
with pin-ended supports. The length of the shaft between supports was 
taken as 40.5 inches. From elementary beam theory, K = 48 ene lens thy 
ibs. /in. The calculated value was 65.33 lbs. /in. 


The spring properties of the shaft were also evaluated experimentally. 
Tests were made with the shaft on knife-edged end supports and with the 
shaft supported in the bearings. Any change in response would indicate 
a resisting moment could be exerted despite the self-aligning feature of 
the bearings. Tests were made for the unloaded shaft and with the rotor 
in position at midspan. With the shaft on Supports placed 40.5 inches 
apart loading and unloading was conducted in the same manner described 
for the springs. Results are listed in Tables IV and V for the knife- 
edge and bearing support end conditions respectively. Figure 6 shows 
the experimental results graphically. No significant change resulted 
from the change in end support conditions. From the slope of the curve 
in Fig. 6 the shaft stiffness was calculated as 66.5 lbs./in. Although 
the test with the bearing showed some hysteresis loss the effect was 
small, and was therefore disregarded. The experimental value, in close 
agreement with the theoretical value, was used in all subsequent cal- 


culations. 


She 





Table IV 


STATIC SHAFT SPRING CALIBRATION - KNIFE EDGE SUPPORTS 











Mor 
Wt. (lbs) Deflection(in. ) МЕ. (1Ъ$.) Deflection(in. ) 

0.0 0.0 и +0.0 0.0 
1.0 0.014 "+1.0 0.016 
2.0 0.029 +2.0 0.032 
220 0.042 Ea 0.046 
4.0 0.057 +4.0 0.064 
5.0 0.072 +5.0 0.080 
6.0 0.088 +6.0 0.090 
7-0 0.103 +7.0 0.107 
8.0 0.116 +6.0 0.093 
9.0 0.132 +5.0 0.079 
10.0 0.146 +4.0 0.064 
11.0 0.160 +3.0 0.048 
12.0 0.178 +2.0 0.032 
11.0 0.164 +1.0 0.016 
10.0 0.147 W +0.0 0.003 
9.0 0.12 B 

8.0 0.117 

720 0.102 

6.0 0.086 

5.0 0.072 

4:0 0.056 

3.0 0.042 

2-0 0.029 

1.0 0.014 

0.0 0.0 

*W = 14.5 lbs. 
r 
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Table V 


STATIC SHAFT SPRING CALIBRATION ~ BEARING SUPPORTS 


*W 


rotor 
Wt (1bs.) Deflection(in.) МЕ. (1Ъ$.) Deflection(in. ) 

0.0 0.0 и +0.0 0.0 
1.0 0.015 о 0.019 
2.0 0.028 +2.0 0.041 
320 0.041 +3.0 0.050 
4.0 0.061 44.0 0.070 
5.0 0.076 45.0 0.084 
6.0 0.091 46.0 0.103 
210 0.108 +7.0 0.115 
8.0 0.125 +6.0 0.103 
9.0 0.145 +5.0 0.090 
10.0 0.157 44.0 0.074 
11.0 0172 +3.0 0.060 
1250 0.190 +2.0 0.046 
11.0 0.181 NE 0.030 
10.0 0.164 W +0.0 0.014 
9.0 0.149 

8.0 0.15 

7.0 0.121 

6.0 0.103 

5.0 0.088 

4.0 0.076 

3.0 0.055 

2.0 0.042 

1.0 0.029 

0.0 0.016 


*W = 14.5 lbs. 
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3.2.4 Rotor and Lip Weight 
Details of the rotor are shown in Appendix A and listed in Table Т. 


Figure 7 shows the relative size of these parts. 


The rotor was balanced statically. Dynamic balancing was considered, 
but it was felt that the disk, only 1.5 inches thick, would produce 
only small effects due to possible variations in properties between its 


parallel planes. 


The rim weight, which was introduced after static balancing, was 
employed for two purposes. Knowledge of its size permitted calculation 
of rotor eccentricity, defined as the shift in center of gravity of the 
rotor due to the added weight. This value is essential to accurate 
calculation of exciting forces developed during rotation. Secondly, 
considering the aforementioned static balancing, the imbalance due to 
the lip weight would overshadow effects of any small variations not 


removed in the balancing. 


The eccentricity was calculated as follows: 


W = basic rotor weight, lbs. 
w = added rim weight, 0.04 lbs @ 2.75 in. from axis 
x location of displaced c.g. from polar axis, in. 


li 


X 


Moment about Xx = W X - w(2.75 - x) = 0 


c 0.04 x 
3 14.44 + 


t 


TE 


0% = 0.0076 in. = € 





© 


3.2.5 Control System 
Details of the components of the control system are shown in 


Appendix A. The general arrangement is shown in Fig. 2. 


The item referred to as the control valve in Fig. 2 was not designed 
specifically for this use. It is identical to the needle valves used as 
metering devices for the working fluid in the dashpots. Large adjustments 
in supply of compressed air to the drive jet were made by means of a 


coarse valve in the line upstream of the control valve. 


3.2.6 Dashpots 
Detailed drawings of the dashpots are shown in Appendix A. 
Smith, | Ди has written, 


The assumption as to the viscous law of damping is not fulfilled 
in practice and analytical predictions as to the effect of damping 
are, to some degree, only qualitative. 
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The dashpots required extensive testing to determine their char- 
acteristics. It was desired to determine three main factors: 
(1) the degree to which the dashpots were linear, (2) the dependence 
of the magnitude of damping on the direct communication port openings 
(See Fig. 8), and (3) the range of damping "constants" possible 


within limits of needle valve settings. 


Any attempt to evaluate the dashpots by analysis of the internal 
flow would require numerous assumptions. An analytical description of 
the non-steady flows through the bypass tubing and direct communication 


ports would have resulted in a doubtful approximation at best. 


A method based on direct measurement of drive force, acceleration 


and velocity response was used. 


3.2.6.1 Instrumentation and Test of Dashpots 

A dashpot was mounted on a fixed platform above a Ling Electronics 
Model 219 Shaker as shown in Fig. 9. А drive spindle of 6061-T6 aluminum, 
{8}, was made to accept the male thread of the dashpot piston shaft. The 
spindle, except for the end portions, was square in cross section. The 
width of the sides of the cross section was machined to 0.250 + 0.001 
inches. Two SR-4, type A-7 strain gages, {9}, were mounted axially on 
opposite sides of the spindle. The gages were connected in series to 
increase sensitivity to strain, and at the midpoint of the spindle to 
avoid effects from the non-uniform ends. The other end of the spindle 


was threaded into an adaptor on the drive table of the shaker. 


Two accelerometers were mounted on the adaptor. Figures 10 and 11 
illustrate the details of the arrangement. A Statham Model 5A5-100-120 
accelerometer was used to meaSure accelerations. Its output, as well 
as that from the strain gages, was led to the Honeywell Model 130-2C 
Carrier Amplifier shown as item 7 in Fig. 12. Both amplified signals 
were recorded semi-permanently on the output paper of the Heiland 
Model 906B Oscillograph Visicorder, item 6, Fig. 12. An electronic 
timer imposed time reference lines on the simultaneous trace of accel- 
eration and strain. The output of the second accelerometer, a Consol- 


idated Engineering Corporation Model 4-102-A Vibration Pickup, measured 
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velocity and peak-to-peak amplitude of vibrations. Readings of average 
velocity and displacement were indicated on a Consolidated Electro- 
dynamics Corporation Type 1-117 Vibration meter. The velocity, amplified 
by the vibration meter, was led to a Hewlitt-Packard Model 120A 
Oscilloscope. A typical trace of the variation of velocity with time 

is shown in Fig. 13. Operating frequency was read directly from the 


frequency control dial on the shaker control console. 


The strain trace on the Visicorder was calibrated in accordance 
with operating instructions for the system. The dashpot was charged 
with H-515 OHA petroleum base hydraulic fluid, {10}. Several combin- 
ations of port openings and bypass valve settings were observed. 
Operation with hydraulic fluid in the dashpots was conducted in some 


cases and air was used in other instances. 


lt was believed that even the minimum possible damping on the part 
of the dashpots might prove excessive in relation to other characteristics 
of the system. It was decided to attempt to approximate the case of 


zero damping. 


The dashpot configuration which would produce minimum damping 
obviously required the bypass line and communication ports to be unre- 
stricted. The only question lay in the choice between air and hydraulic 
oil as the working fluid. Although the choice of air would seem the 
more logical, there was a higher likelihood of large frictional forces 
associated with its use. It is not known whether an air loaded dashpot 
would, in fact, produce less damping than that which was measured, but 
the hydraulic fluid was selected for use in the test based on its lubri- 


cating properties. 


The plugs were removed from the communication ports and the bypass 
valve was set fully open. Average velocity for the test was 3.0 in./sec. 


Peak-to-peak displacement was 0.28 in. at the test frequency of 5 cps. 


Traces of acceleration and strain (force) were scaled up four-to-one 
for readability. Figure 14 shows only the basic variations of force 
and accelerations as recorded on the trace. Superimposed on the actual 


trace was a high frequency "hash". This was present during the warm up 
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of the shaker prior to actual tests. Since these signals were obviously 
not associated with response of the dashpot they were disregarded in the 
interpretation of the test traces. By successive integrations of the 
acceleration trace, velocity and displacement traces were produced. 

These waveforms are shown in Fig. 14 a-c. The associated force trace 

is shown in Fig. 14 d. Integration of the acceleration wave produced a 
parabolic velocity trace. The velocity trace was scaled from the known 
value of the average velocity, two-thirds of the maximum velocity repre- 
senting the averageheight of the parabolic trace. The scale of the 
acceleration trace was found from the fact that the area under the accel- 


eration curve was equal to the known maximum velocity. 


According to design specifications the shaker should have produced 
a sinusoidal output with maximum of 3 percent distortion at 5 cps. The 
derived displacement-time relationship was not a sine wave. Whether this 
was due to instrumentation error or to shaker operation with distortion 


greater than that cited in the specifications is not known. 


With the assumption of sinusoidal motion for the system depicted in Fig. 15 


the following relationships apply: 








F. = inertia force =m x : Fa = damping force = сх 
Fr = total force = Fi T Ey 
. . .. 2 | 
X = A sin Ut 5 x = WA cos Wt ; x = - 0 À sin E 
NE SURE. pus IER X 
Sinot = д = | 29) созба ош | (3) 
Ж Er) еы 
F+ = mw A UR ) EE t Са — AN ELLIPS 


Баден с xX , a straight line whose slope is c. 


2 
FI = -mW x , a straight line whose slope is - т we, 
А-/ x Ис ° ВИ 1 › АМ ELLIPSE 
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Figure 15 shows these relationships graphically. The ellipse in the 
velocity plot has a minor semi-axis proportional to the inertia or 
mass of the system. The minor semi-axis of the displacement plot is 
proportional to the force of viscous damping in the system. Repre- 
sentative plots of experimental data are included for purposes of 


comparison. 


Consistent with the nature of the forces involved the inertia force 


contributes nothing to the area of the displacement plot. 


Test data plotted in Fig. 16 indicate the variation of total force 
with velocity. The shape of the enclosed area differs markedly from the 
ideal form of an ellipse. The distortion of the figure emphasizes the 
remarks cited at the beginning of this section. The non-linearity is 
strongest in the region of negative velocity. Deviation from the ideal 
case is also apparent in the recorded force trace, Fig. 14 d. This may 
have been caused by coverage of the bypass ports at the extreme point of 
the piston stroke, by large friction forces present due to momentary 
loss of lubrication between the dashpot cylinder and the piston surfaces, 
or by lag in the mechanical linkage between the drive table and the 
piston. The vertical offset of the center of the closed figure from the 
origin indicates that the force is dependent to some degree on the 
absolute direction of motion. No logical reason for such a dependence 


is apparent. 


For a linear system the slope of the mean line of the F, - velocity 
plot is equal to the damping constant. Since no dependable mean line 
is apparent, the slope of a line drawn through the endpoints of the figure 
is considered only as a first estimate. The slope so defined was cal- 


culated as 1.88 lbs./(in./sec.) 


In an attempt to analyze the Fr ~ velocity trace further, the product 
of mass and acceleration was subtracted point-by-point from the total 
force ordinate. Figure 17 shows the result as a "Reduced Force"-velocity 
plot. A further refinement of this figure could be made by an allowance 
for the mass of fluid being accelerated during the operation of the 
dashpot. Regardless of the distortions in Figs. 16 and 17 it can be said 


that, essentially, the dashpot response has been established. 
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The force-displacement plot from test data is shown in Fig. 18. 
The closed curve has some resemblance to an ellipse. Assuming that 
an equivalent damping system will absorb or dissipate equal energy per 
cycle as does the physical system, an equal area ellipse was constructed. 
One semi-axis was defined by the known displacement. From the formula 
for the area of an ellipse, Area = Tm xa x b, the regular ellipse, show 
dashed in the figure, was drawn and superimposed. From the scale of the 
test figure the semi-axis proportional to cWA was measured. Dividing by 


the known value of WA the value of the damping constant was found. 


5 


толе лкы ss 1.77 1bs. /(in./sec.) 


с = 


This value, in fair agreement with that calculated from Fig. 16, is 


probably the better of the two. 


3.3 Critical Speed Test Program 


Tests of the rotating system were conducted with the bearings locked 
and with the bearings under the restraint of the dashpot-spring combination. 
The dashpots for the latter test were set in the same configuration used 
in their calibration, i.e., direct communication ports and bypass control 


valve fully open. The intermediate spring was employed with the dashpot. 


3.3.1 Instrumentation 

The electrical probe displacement sensor shown in Fig. 19 was used 
to measure deflections. A conducting probe in the form of a small 
spring was fixed to the end of the vernier shank. The vernier micrometer 
was mounted rigidly in a support. A 22.5 volt dry cell battery provided 
potential to the probe through the conductor shown. This voltage was 
impressed in parallel across the vertical deflection circuit of the 120A 
Oscilloscope. At each instant of contact between the probe and the 
rotor the scope circuit was grounded and a vertical spike appeared on the 
scope. Zero readings of the vernier were made prior to each test. The 
vernier, readable to 0.001 inches, was then backed off. For each speed 
at which measurement of vibration amplitude was desired the vernier was 
advanced until first contact was indicated on the scope. Some arcing 
at the tip of the probe indicated small errors, estimated at + 0.005 in., 


were inherent to the method. 
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An RCA Mode} 6318 Strobotac was used for speed indications. The 
low scale on this instrument, usable for speeds from 60 to 3700 RPM, was 


readabie to + 2 RPM. 


3.3.2 Results 
At no time during tests with the bearings supported by the inter- 
mediate spring was there any perceptible or measurable motion of the 


bearings. 


Results are listed in Table VI, Experimental Results. Results are 
shown graphically in Figs. 31 and 39 in conjunction with the comparison 


of theoretical and experimental results. 


3.4 Comments on the Physical System 
The system could not be tested without the dashpot linkage, and 


consequently, the dashpots engaged. This resulted from the inability of 
the spring-spring holder combination to resist loads normal to its 

axis. The stabilization provided by the dashpot linkages was desirable 
in one sense and indicative of a possible source of trouble in another. 
Forces along the rotor shaft axis were carried through the dashpot 
linkages to the dashpot pistons. An aggravated case of this could have 
resulted in the pistons becoming cocked in the dashpot cylinders. Such 
cocking could have caused deviations from linearity even greater than 


those observed in the controlled testing of the dashpots. 


The light spring was not capable of Supporting the combined weight 
of its associated bearing and the rotor mounted on the shaft. These 
springs "bottomed" in their holders when the fourteen pound rotor was 
mounted. A lighter rotor, capable of being supported by the light 
springs, was mounted. Spring action was linear with the rotor at rest 
but, movement of the bearings in vibration closed the springs coil-to- 


Coll. 


The drive disk whipped noticeably during strong vibration of the 
system. Air was spilled from the cusps in the disk and a non-steady 
drive resulted. The disk could be reworked to form channels as in a 
Pelton type turbine; this would reduce spilling. This same modification 
would eliminate the small axial component associated with the present 
drive system. 


On 
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Table VI 


EXPERIMENTAL RESULTS 


Bearings fixed Hydraulic fluid in dashpots, 
bypass and ports open 

RPM Deflection (in.) RPM Deflection (in.) 
250 0.044 300 0.001 
260 0.045 330 0.000 
270 0.046 340 0.005 
280 0.047 335 0-029 
290 0.048 370 0.040 
300 0.050 Зе 0.080 
310 0.062 380 0.104 
320 0.063 383 0. 116 
330 0.066 388 0. 109 
340 0.067 390 9113 
350 0.090 E 0.090 
360 0.099 400 0.075 
370 0.238 450 0.058 
1/2 0.670 

375 0.969 

380 1.245 

385 0982 

390 0.07? 

395 0.051 

400 9:017 

450 0 005 

500 0. 003 


- 39. 





и. 





w. p we 4 





ñ 


ri 717: 


Strong transient vibrations were noted during accelerations of the 
rotor. The speed control system is not conducive to setting of precise 
speeds. Low speeds require small inputs from the jet. Acceleration forces 
must be kept very close to those necessary for steady drive at the speed 
desired or the speed will build rapidly and overshoot. For low speeds 
this means allowance for long periods of speed buildup. Fifteen minutes 


was not excessive for speed stabilization during tests. 


Fluctuations in the supply air pressure caused speed to drift even 
after long periods allowed for stabilization. In consideration of these 
factors measured speeds can be interpreted as accurate only to within 


+ 5 RPM. 


Error in the calculation of the dashpot damping constant may have 
stemmed from either or both of the following: (1) Low amplification of 
the spindle strain trace such that its maximum value reached only twenty 
percent of the lowest mark on the calibration scale of the recorder. 

(2) Reading and scaling of recorded traces. These sources of possible 
error could be eliminated or minimized by, (1) use of more sensitive 
Strain gages with larger output per unit strain, and, (2) recording 
velocity waveform in addition to simultaneous acceleration and drive 
Spindle strain. This would require the use of a strain gage type 
vibration pickup; the Consolidated instrument used in the tests described 
was a self-exciting, piezo-electric type whose output could not be 


matched to the recorder galvanometer. 


Other configurations of bearing support, in particular that with the 
light spring mounted in the horizontal plane, were observed. Bearing 
motion was pronounced and the spring appeared to be operating without 
bottoming. This constituted a system with non-symmetrical bearing 
stiffness. Although not within the scope of this study, this configuration 


is of considerable interest any may form a logical extension to this work. 
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4. THEORETICAL ANALYSIS 


A model representative of the rotating system is shown in Fig. 20. 
Component names and general nomenclature associated with the physical 


system were retained wherever possible. 


Values of the physical properties of the system used in this 
analysis differ somewhat from those listed in Table I. Only three 
masses are defined in the theoretical model; the mass of the rotor, 
m, > and the masses of the left and right-hand bearing complexes, m, 
and та) respectively. 

4.1 Equivalent Mass Lumping 

Any physical system can be said to have infinite degrees of freedom 
in that it is composed of an infinite number of point masses. In some 
cases less significant masses are ignored completely. An example is the 
elementary spring-mass system. The center of gravity of the spring is 
obviously displaced during oscillations and it therefore constitutes a 
mass in motion. Since the principal mass of the system is very often 
several orders of magnitude greater than that of the spring and the 
displacement of the spring mass is so small its contribution to the 
system is insignificant. For the system under consideration here, an 
attempt is made to formulate a logical expression to reflect the con- 
tribution of the motion of the shaft to the total energy of the system. 
The kinetic energy of the actual shaft in vibration will be referred to 
as the shaft energy. The shaft energy will vary for each mode of 
vibration. Figure 21 shows, in a representative manner, the three modes 
of vibration for the three mass system. It is easily seen that no 
single factor will suffice in apportioning the mass of the shaft among 
the masses in the model. The first mode shown requires that the con- 
tribution be made solely to the rotor. The other modes indicate that 
the bearing masses should be increased by some factor but that the rotor 
mass need not be modified. If the elastic line of the shaft is assumed 
to be parabolic, eight-fifteenths of the shaft mass should be added to 
that of the rotor for similarity. If the other modes were calculated 


on this basis they would be in error. Ideally an assumption should be 
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made as to the apportionment for each mode, the mode shapes calculated 
and the respective factors revised. Due to the relatively small mass 
associated with the actual shaft an iteration of this kind was not felt 
to be necessary. A factor of one-third was taken; the mass of each 
bearing and that of the rotor was increased by one-third the mass of the 
adjacent shaft section(s). The mass of the overhung drive disk, being so 
physically close to the left-hand bearing, was lumped with that bearing 


complex. 


Table VII, Equivalent Lumped Masses, lists the calculated values 


of the three masses in the theoretical model. 


4.2 Mass Dynamics 


Free body diagrams of the three masses with forces acting during 
rotation are shown in Fig. 22. The differential equations of motion 
based on the free body diagrams reduce to a system of six simultaneous 
algebraic equations. The development is shown on the pages following 


Bio. 22. 


Had a solution in the form X, = Ret been assumed, the system 
could have been described in three equations with complex coefficients. 
Computer programming of complex expressions requires a somewhat more 
demanding technique, and, for this reason, the expressions with real 


quantities throughout were selected. 


The six equations describe the theoretical model of the rotating 
system in damped, forced vibration. As shown in the last page of the 
development, the system of general equations reduce to a specific case 


for the rotor at the middle of the shaft span. 


4.2.1 Program CRTSPD 

The six simultaneous algebraic equations are the basis of the 
computer program CRTSPD (Critical Speed). Fortran ,{11}, coding was 
used. A program flow chart indicates the logic and progress of the 
computations in the computer. A sample program is included in 


Appendix B. 
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Rotor mass: 


m, = ue + = (.0544/386) (20.25) = 0.394 





Left-hand bearing complex mass: 


m, = == + ج‎ (.0544/386) (20.25) + 2.035 = 0.0270 





Right-hand bearing complex mass: 


m, = 7.88 + 1 (.0544/386) (20.25) = 0.0220 


Table VII 


EQUIVALENT LUMPED MASSES 


Significant Mass - W(lbs.) M(m.u.) 

L.H. Bearing complex 10.434 020270 mq 
Rotor 152218 0.0394 т) 
R.H. Bearing complex 8.501 0.0220 m, 


Mass Ratios: 


0.0270/0.394 = 0.6858 


m, /m 
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m./m, 0.0220/0.394 = 0.5588 
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Figure 23 


Flow Chart PAROUCRAM CuT5PD 


STARTE 
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The dimension statement apportions storage space for the arrays 
representing the augmented matrix of coefficients of the algebraic 
equations of motion and the values of their solution, X(6). Column 
headings for the printout are specified. The constants of the physical 
system, Cl-C6, represent We ,“س‎ шц“, E Mz, and e, respectively. 

CA is the coded form of the damping constant at the left-hand bearing. 

CE is the corresponding damping constant for the right-hand bearing. 

These constants were set equal to one another for all calculations. The 
value so assigned is hereafter referred to as the symmetrical damping 
constant. The DO 10 statement defines the range of RPM. Omega is then 
defined in terms of RPM. Subroutine CALCA, to which all factors necessary 
for the calculation of the elements of the augmented matrix have been 

made available through use of the COMMON statement, then calculates each 
element according to the specified value of RPM. The resulting numerical 
array, A, common to Subroutine JORDAN2, is solved as a set of simultaneous 


equations. The values of the solution are returned to the main program 


and amplitudes of mass motion and associated phase angles are calculated. 


Calculated system response was initially investigated over a broad 
range of speeds. Results indicated that the critical speeds of the 
system lay between 200 and 900 RPM. Calculations for unit steps of RPM 
were selected so that variations of response with speed would be well 


defined in the critical regions. 


Values initially chosen for the damping constant were somewhat 
arbitrary. Calculations for very small and very large values of the 
damping constant were made. Investigation of system response for all 
possible combinations of damping constants at the two bearings was obvi- 
ously impractical. The case in which the two constants would be kept 
equal seemed to be a suitable starting point. Hereafter the term 
"symmetrical damping constant" or merely "c" will be used to describe 


the stated relationship. 


From the preliminary computer calculations it was noted that the 
most profitable region for close investigation was defined by small 


values of c, i.e., between 0.0 and 3.0 lbs/(in./sec.). 
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Computer runs were made for no damping, (bearings supported by the 
springs alone), and for damping values up to 3.0 1bs./(in./sec.) in steps 
of 0.5 lbs./(in./sec.). The case, supposedly analogous to the experimental 
test for damping of 1.77 lbs.(in./sec.), was also run on the computer. A 
value of exceptionally large damping, 10,000 1bs./(in./sec.), was assigned 
to produce calculated response corresponding to the fixed-bearing exper- 
imental test. Results of the computer calculations described are shown 
graphically in Figs. 24-41. The three modes of vibration associated 


with no damping are shown in Fig. 24. 


4.2.2 Critical and Optimum Damping 

Ordinarily the concept of critical damping is reserved for systems 
with few degrees of freedom in free vibration. In the simple, single 
degree of freedom, spring-mass-dashpot system as treated for example by 


Burton, а the differential equation of motion is developed as follows: 


.. ә des 
MX == CX-*Rx ASSUMING 4-28 57+ Е б +: =о 
Mm 
5 < -C Ł VC*- 4m ё, 
2 77u 

Бов сто, FOR CRITICAL DAMPING, 

С? = 4m k. C= 2Vm #. 

^ 25 


X= Pe 


The finite value of c for which the term in the radical reduces to zero 
is referred to as the critical damping constant. From the resulting 
expression for the response of the system at critical damping it can be 
seen that the system does not oscillate; any initial displacement from 
equilibrium is followed by a direct return of the mass to the equilibrium 
position. Larger values of damping result in a similar response but the 
time to return to the equilibrium position is increased. For infinite 
damping the time for return is infinite, or, in other words, the mass 
never returns to its initial position. The pertinent part of the defin- 
ition of critical damping is that it specifies the minimum value of 
damping for which the mass does not oscillate in response to a disturbance. 


For physical systems damping is always greater than or equal to zero. 
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The single degree of freedom system is described by a quadratic 
in 5. If this treatment is extended to a system with two degrees of 
freedom with both masses subject to damping a quartic expression in 6 
results. In this case two damping constants govern the nature of the 
imaginary parts of the two pairs of complex roots. Generally then, for 
a system of n masses there will be an expression of degree 2n in 9 and 
n damping constants will govern the nature of the n conjugate pairs of 


complex roots. 


Applying this treatment to the system under investigation results 


in a sixth degree expression in 6. 


= а JLT = pt, + Us 2 я Аз = V=00/+0e/4 Cu, “= wi = De 


p. À Е (ar e EE e N)S (2 V Lus) e (V Eu ENS Es — 


о 206 
2, /6 E 


Certain properties of the roots can be predicted from physical consid- 
erations alone. Since only two of the masses are subject to damping 
there will be only two governing damping constants. One mode of oscil- 
lation will be present no matter what values of damping are assigned. 
For certain values of damping the imaginary parts of two of the three 


complex pairs of roots will vanish. 


An extensive analysis of the sixth degree expression is beyond the 
scope of this study. The computer was used in a limited analysis. Results 
indicated that, for the masses and stiffnesses involved in the present 
tests, the imaginary parts of two of the pairs of roots would vanish for 
damping constants of about 3.0 and 5.0 lbs./(in./sec.) respectively. For 
these values the system, in free vibration, would exhibit a response some- 


what analogous to a less complex system with critical damping. 


Figures 24-41 show the calculated response of the system over a range 


of small damping constants and for one extremely large value. 
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Figure 24 shows the modes of the three critical speeds for the 
undamped case. The calculated amplitudes of the masses at these three 
frequencies are not infinite because the exact values at which they occur 


were not assigned in the speeds selected for calculation. 


Ав c is increased from 0.0 to 2.5 lbs./(in./sec.) the amplitude response 
of the masses at the first critical speed is seen to decrease. See 
Fig. 41. Responses at speeds higher than the first critical speed become 
less and less pronounced as damping is increased. Increasing the damping 
has a relatively greater effect on the amplitudes at critical speeds 
higher than the fundamental; once the damping exceeds 2.5 lbs./(in./sec.) 


there are virtually no critical speeds higher than the fundamental. 


Comparison of Figures 35 and 37 shows a reversal in one trend of the 
response for increasing damping. Although the amplitude response of the 
bearings continues to decrease with higher damping, the rotor response 
is larger at 3.0 than at 2.5 lbs./(in./sec.). For the strongly damped 
case shown in Fig. 39 the rotor response has increased greatly over that 


for the lightly damped cases. 


Figure 41 shows the amplitude response of the masses as a function 
of the damping constant. It can be seen that the value of 2.5 lbs./ 
(in./sec.) constitutes that value of c for which the amplitude response 
of the system is a minimum. From this it is apparent that, under these 
conditions, maximum energy is being dissipated from the system. This 


case is defined as that of optimum damping. 


4.2.3 Discussion of Theoretical Results 

Figures 40 and 41 show the major effects of damping on the system. 
Figure 40 shows that the critical speed depends directly on the magnitude 
of damping. Similarly, Fig. 41 shows the dependence of amplitude response 


on damping. 


The limiting case of damping on critical speed is that for which 
damping is infinite, or, in other words, fixed supports. The system is 
then reduced to one of the single mass of the rotor. The calculated value 


of this critical speed is 392 RPM. 
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In order to complete the description of the system at each value 
of damping, phase angle diagrams are included with each figure showing 
amplitude response. Somewhat contrary to the behavior of the spring- 
mass-dashpot system, the system considered shows a marked dependence of 
critical speed on damping, and the amplitude response is not always in 


inverse proportion to the damping. 


4.3 Comments on the Theoretical Analysis 


In general, the theoretical results appear to be valid and consistent 
with established theory. In confirmation of this, a simple, if limited, 


check on the method of calculation is presented. 


For the case of fixed supports, Prohl, | . 


Ni = first critical speed, RPM 


k 


u 


shaft spring constant, lbs./in. 


W = rotor weight, lbs. 


N 


3 
| = 187.7 (k/W)* RPM 


5 
187.7 (66.5/15.21)* = 392.3 RPM 


i 


The calculated critical speed for this configuration was 392 RPM. 


E 





5. COMPARISON OF THEORETICAL AND EXPERIMENTAL RESULTS 


Figure 31 shows the calculated and test amplitude responses for the 


case of damping of 1.77 lbs./(in./sec.). 


The observed critical speed was 383 RPM. From Fig. 40, First 
Critical Speed vs. c, the critical speed corresponding to c = 1.77 lbs./ 
(in./sec.) is 346, a discrepancy of 37 RPM. Also from Fig. 40, a critical 
speed of 383 RPM corresponds to a damping factor of 6.5 lbs./(in./sec.). 
From analysis of the observed critical speed there is the indication that 
the damping during the test was greater than that of 1.77 lbs./(in./sec.) 


thought to have been in effect. 


Amplitude indications from the test were analyzed in conjunction with 
Fig. 41, Mass Amplitude at First Critical Speed vs. c. From the figure 
it can be seen that an amplitude of 0.116 corresponds with either of 
two damping constants, 0.047 or 18.0 1bs./(in./sec.). The lower value is 
felt to be totally improbable. The higher value, although excessive in 
itself, tends to confirm the indication noted in comparison of critical 
speeds that the damping was larger than 1.77 lbs./(in./sec.). Amplitude 
measurement error, although possibly greater at critical speeds due to 
transients from non-steady drive, could not account for the large dis- 


crepancy. 


The apparent discrepancy in dashpot damping could stem from any of 
several sources. The dashpots may have been cocked during the run 
and produced a strong non-linear response. The dashpots were mounted 
vertically for their calibration. Their operation in the horizontal 
plane during the critical speed test may have produced the apparent 


change in their characteristics. 


Figure 39 shows the calculated and experimental response for the 
fixed-bearing case. The calculated critical speed was 392 while the 
experimental value was 380. Since amplitudes are theoretically infinite 
for the fixed-bearing case, comparison is not applicable as in the case 
with lightly damped bearings. Based solely on the comparison of critical 


speeds the agreement between calculated and test results is fair. 
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6. CONCLUSIONS 


The system of a rotor mounted on a shaft with damped and flexibly 
Supported end bearings shows its amplitude response to be dependent on 
bearing support impedance. Rotor amplitudes decreased as damping in- 
creased to the optimum value of 2.5 1bs./(in./sec.) and increased for 
larger values of damping. There is a limit to which rotor amplitude 
at critical speeds can be influenced by damping at the bearings and one 


critical speed will persist in spite of damping applied. 


The calibrated damping characteristics of the dashpots were either 
in error or dashpot damping was effected by the mechanical arrangement 


of the dashpots on the critical speed rig. 


Accuracy of test measurements could be improved. 


ль 





7. RECOMMENDATIONS 


The critical speed test rig offers promise for use in further 
vibration analyses if improvements in its damping and control charac- 


teristics can be effected. 


It is recommended that the physical system be improved by 
a. Modification of the drive system to eliminate thrust along the shaft 
axis, 
b. Design of a support for spring-spring holder assembly to permit testing 
without dashpots, 
c. Provision of springs with stiffnesses between those of the springs 
tested, and 


d. Use of better speed and amplitude indicators. 


It is also recommended that a further investigation of the damping 


range and characteristics of the dashpots be conducted. 
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